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We present a non-perturbative study of the massive Schwinger model. We use a Hamiltonian 
approach, based on a momentum lattice corresponding to a fast moving reference frame, and equal time quanti- 
zation. 



1. INTRODUCTION 

The most important non-perturbative approch 
to study QCD is lattice gauge theory. How- 
ever, there are some observables, where com- 
putational progress in lattice gauge theory has 
been slow. Examples are: Finite density ther- 
modynamics (quark-gluon plasma), exited states 
of the mass spectrum (hadrons and mesons), dy- 
namical scattering calculations of cross sections 
and hadron structure functions, in particular in 
the region of small Q 2 and small xb (10 -2 to 
10 -5 ). A non-perturbative Hamiltonian approach 
may be a viable alternative. The Hamiltonian 
method does have some advantages over the La- 
grangian method. It is relatively simple to ob- 
tain wavefunctions of the hadronic states, and it 
is likely that some physical information on the 
glueballs, e.g., the wavefunction can be derived in 
this formulation more easily . In deep inelas- 
tic lepton-hadron scattering, the success of the 
parton model suggests the physical idea to use a 
fast moving frame also for a computational study 
of those processes. The parton model can be jus- 
tified using the operator product expansion. In 
equal-time quantization, the Breit-frame is the 
most convenient choice of frame in order to in- 
terpret structure functions. In Ref. || the authors 
have proposed such a scheme based on equal-time 
quantization, using a lattice Hamiltonian on a 
momentum lattice corresponding to a fast mov- 
ing frame (Breit-frame). It has been applied to 
the scalar </> 4 theory in 3+1 dimensions. Distribu- 



tion functions and the mass spectrum in the close 
neighbourhood of the critical line of the second 
order phase transition have been computed. Here 
we study a simple gauge theory, namely QED 1+1 

01 , the so-called massive Schwinger model, in this 
framework The purpose of this work is to 
show: (a) The use of a fast moving frame, such 
that v < c, in conjunction with equal-time quan- 
tization works well also in QED\ + i. We obtain 
quite precise results in the ultra-relativistic region 
m/g — » 0. (b) We consider as most important 
new results of this work the non-perturbative re- 
sults of the dependence of the vector and scalar 
mass on the 0-angle. 

2. Method 

Starting from the Lagrangian, we use the axial 
gauge A 3 — to obtain the Hamiltonian 

// - [ dx 3 (^ 3 id 3 ip +iptp) 

,2 rL 



+ 9_ j dx 3 ( ^)_L ( ^). 



(1) 



One introduces a space-time lattice given by spac- 
ing a with N = — lattice nodes. Via discrete 
Fourier transformation one goes over to a mo- 
mentum lattice, with cut-off A = — and rcsolu- 
tion Ap — Motivated by the parton picture 
we make the assumption that left-moving parti- 
cles are not dynamically important, if physical 
particles are considered from a reference frame 
characterized by a velocity v = which is close 
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Figure 1. The dimensionless mass of the vector 
boson versus m/g in the chiral region. 



Figure 2. Dependence of the dimensionless mass 
of vector and scalar particle on the (9-angle after 
subtraction of t9-vacuum energy. 



to the velocity of light. We thus consider a mo- 
mentum lattice where < p°,p 3 < P. In order 
to minimize the number of virtual particle pairs 
created from the vacuum, we choose a small lat- 
tice size. The reason for this is that the number 
of vacuum pairs is roughly proportional to a vac- 
uum density times the lattice size. On the other 
hand, a fast-moving physical particle is Lorentz 
contracted; thus it fits into a small lattice vol- 
ume (when compared to the rest- frame). For 
the purpose of computing the mass spectrum, we 
need to determine the vacuum energy. Because 
the vacuum has the quantum number P = 0, 
the vacuum energy (and only this) is computed 
in the rest frame. Because the model is super- 
renormalizable, one can perform the continuum 
limit a — > 0. The only renormalization necessary 
is the subtraction of the vacuum energy. On a 
space-time lattice, one has to satisfy a physical 
condition (scaling window) a < £a < L, where 
£ is the correlation length in dimensionless units, 
related to the physical mass of the ground state 
by M = j^. In a strongly relativistic system, 
M <C P, thus when a — > the scaling window is 
replaced by X/P < 1/M < L. For more details 
compare with Ref . || . 

3. Numerical results 



3.1. Mass spectrum 

We diagonalize the Hamiltonian in a sector 
with momentum P — to obtain the vacuum 
energy E vac . Then we diagonalize the Hamilto- 
nian in a sector P ^ corresponding to a rel- 
ativistic velocity. This yields an energy spec- 
trum E' n . The physical energies are obtained 
from E n = E' n — E vac . The mass spectrum is 

then given by M n = Je% - P 2 . The low ly- 
ing states of the massive Schwinger model are 
a vector state and next a scalar state ||. The 
mass of the vector boson in the chiral region is 
shown in Fig. [l| The vector particle is almost 
entirely a fermion-antifermion (qq) bound state. 
We find good agreement with chiral perturbation 
theory [Q and with finite lattice results by Hamer 
et al.§. 

3.2. Dependence on (9-angle 

The massive Schwinger model has 9- vacua and 
one can study its fraction. The wavefunction 
is invariant under local gauge transformations. 
The t9-angle characterizes the behavior of the 
wavefunction under global gauge transformations 
V[A] -» e in9 ^[A], where n = 0, ±1, ±2, • • •. The 
mass spectrum of low- lying states as a function of 
the #-angle is shown in Fig. 0. The results are in 
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Distribution Function of the Vector Boson 

N=192 




Figure 3. Distribution function of the vector bo- 
son for different m/g in the chiral region. 



agreement with first oder chiral perturbation the- 
ory in the ultra-relativistic regime (m/g < 0.04 
for TV = 6 and m/g < 0.01 for N = 24). The 
Feynman-Hellmann theorem relates the fermion 
condensate < ipip > to the derivative of the vac- 
uum energy. Similarly holds for the vector state 



< tpip >- 



M<°) d 
2tt dm 



M, 



(2) 



Extracting the slope 8M v /dm from our data we 
obtain < -tpip > /g = 0.16cos((9) while the exact 
solution of the massless Schwinger model gives a 



factor 



27T\/7r 



0.1599 on the r.h.s. 



3.3. Distribution functions 

In the Hamiltonian approach it is easy to com- 
pute the wavefunction of a low-lying state. ^From 
the wavefunction one can obtain information on 
its parton structure, i.e., the number of partons 
and their momentum distribution. The distribu- 
tion function of the vector boson is given by 



f(x B ) =< *„(P) 



1 



b\b k 



d\d k 



*v(P) >,(3) 



where xb = k/P is the fraction of momentum of 
the vector boson carried by the parton, i.e., In 
Fig. U we display lattice results of the distribu- 
tion function for m/g = to 0.28. In the mass- 
less limit the distribution function has the shape 



of a box. Our results are compared with varia- 
tional calculations using the infinite-momentum 
frame by Bergknoff || and also with those by Mo 
and Perry || using the light-cone. The most sen- 
sitive region is the ultra-relativistic region. We 
find agreement in shape with Mo and Perry's re- 
sults and very good agreement with Bergknoff's 
results. 

4. Summary 

In this work we have applied a Hamiltonian lat- 
tice approach in a fast moving reference frame 
to study the massive Schwinger model. We find 
that the method works well for the computation 
of the low-lying mass spectrum and distribution 
functions. It works also in the presence of the fr- 
action. Here we have investigated only 9 = 0, ir. 
More 0-angles can be studied, e.g. by adding a 
small number of negative momentum states to the 
basis. 
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